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Abstract 



We discuss the generation of the /U-term in a class of supersymmetric models charac- 
terized by a low energy effective superpotential containing a term XSH1H2 with a large 
^c coupling A ~ 2. These models generically predict a lightest Higgs boson well above the 

LEP limit of 114 GeV and have been shown to be compatible with the unification of 
gauge couplings. Here we discuss a specific example where the superpotential has no 
dimensionful parameters and we point out the relation between the generated //-term 
and the mass of the lightest Higgs boson. We discuss the fine-tuning of the model and 
we find that the generation of a phenomenologically viable /i-term fits very well with a 
heavy lightest Higgs boson and a low degree of fine-tuning. 

We discuss experimental constraints from collider direct searches, precision data, 
thermal relic dark matter abundance, and WIMP searches finding that the most natural 
region of the parameter space is still allowed by current experiments. We analyse bounds 
on the masses of the superpartners coming from Naturalness arguments and discuss the 
main signatures of the model for the LHC and future WIMP searches. 



1 Introduction 

At the dawn of the LEP era hopes were on the discovery of electroweak (EW) scale super- 
symmetry, which provided a symmetry principle as solution of the hierarchy problem of the 
Standard Model (SM) and a serious chance for the Higgs boson to be in the reach of that 
machine. The measurements of the gauge couplings performed at LEP turned out to be in 
very good agreement with the idea of a unification at scales Mqut ~ 10 16 GeV, as foreseen 
in the Minimal Supersymmetric Standard Model (MSSM). Unfortunately, the striking picture 
emerging from the study of the gauge couplings did not find a counterpart in the observation 
of either a light Higgs boson or of any of the numerous superpartners. 

However one has to acknowledge that the expectation for SUSY, and in particular for 
a light Higgs boson, to be discovered at LEP was mostly based on strict requirements of 
minimality in model building. The MSSM has in fact a scalar potential entirely fixed by the 
gauge symmetry and this yields the famous result that the lightest Higgs boson of the MSSM 
has to be light. Such lightness of the MSSM Higgs boson can be understood noting that at 
the tree-level the lightest Higgs boson mass can be upper-bounded by 

m z cos 2/3, (1.1) 

and therefore can even be equal to if tan/3 = 1. 

In the MSSM further contributions to the lightest Higgs mass arise at the one loop level 
and can lift the Higgs mass above the experimental lower-bound [1] 

m hsM >114GeV. (1.2) 

However such large radiative corrections require soft masses much larger than the Electro Weak 
Symmetry Breaking (EWSB) scale, which loosens the Naturalness argument to motivate su- 
persymmetry at the TeV scale. 

The necessity of large soft masses in the MSSM motivates the study of supersymmetric 
models that extend the MSSM and that can possibly alleviate the tension with LEP direct 
searches. Indeed a study of extensions of the MSSM with effective operators [2] shows that, 
departing from the minimal model, one can have substantial corrections to the bound on the 



lightest Higgs mass in eq.(l.l). A similar conclusion is drawn in Refs. J3j H] where the issues 
of the Higgs mass and of the unification of gauge couplings have been considered in specific 
renormalizable extensions of the MSSM. 

Refs. El [7j are concrete examples of such models in which the tension with the LEP 
bound is alleviated introducing either extra vectorial fermions or extra gauge structure. In 
these models particular care is taken to preserve the unification of the gauge couplings in 
a manner that is as close as possible to that of the MSSM. However the result is still not 
completely satisfactory because these models either are fine-tuned or mildly affect the maximal 
value of the mass of the Higgs [8]. 

As a matter of fact, relaxing the requirement of strictly perturbative unification in general 
leads to non-minimal models with significant changes for the phenomenology [9j. Indeed, an 
attempt along this line has been made in |10j . where the self-coupling of the Higgs sector 
has a strong coupling phase at some intermediate scale between the EWSB and the GUT 
scale. Handling such a strong coupling phase for the Higgs sector puts one in position to 

1 



substantially increase the mass of the Higgs at the tree-level. The concrete examples in [TO] 
show that this can be done compatibly with the cherished unification of couplings. The idea 
has been further elaborated in [TTJ where it has been discussed a unification-compatible UV 
completion for the so-called Next to Minimal Supersymmetric Standard Model (NMSSM), i.e. 
for a model described by a superpotential of the form 

W = \SH 1 -H 2 + f(S), (1.3) 

where the superfields Hi and H 2 are Higgs doublets and S is a SM singlet. 

With respect to the commonly studied case with perturbative unification, the low energy 
NMSSM of [TTJJ, [TTJ can have larger couplings that become non-perturbative much below the 
gauge coupling unification scale. The extended range of acceptable couplings allows the mass 
of the lightest Higgs boson to reach 200-300 GeV, which leads to dramatic differences in the 
phenomenology of the Higgs sector with respect to the MSSM. 

Indeed, if one considers a model with a superpotential including the term XSHi ■ H 2 , as in 



(1.3), the upper-bound on the lightest Higgs mass is given by 



m% cos 2 2/3 + X 2 v 2 sin 2 2/3 , (1.4) 

which does not vanish for any value of tan /3 and becomes larger as one takes larger values for 
A. 

However the maximal value of A at the EW scale that does not lead to a Landau pole 
below the GUT scale is ~ 0.7 [T2J [T3j, which only modestly affects the upper-bound on the 
Higgs mass, if compared to the MSSM. This motivates the choice to study the regime of the 
NMSSM where the coupling A can be 0(1), namely the regime of the model in which the 
Higgs becomes heavy. Furthermore, the large Higgs mass attainable for A ~ 1 automatically 
reduces the sensitivity of the mass of the Higgs to UV effects, and consequently the model 
turns out to be less fine tuned. This feature adds further motivation to pursue this regime. 



Motivated by this exciting possibility, Ref. |T4] analysed the framework where (1.3) is 
taken as a low energy effective superpotential below a scale of O(10 TeV). Given the great 
importance covered by the coupling A, this framework goes under the name ASUSY. 

In Ref. [H] a very detailed analysis of the Electroweak Precision Tests (EWPT) has been 
performed with the result that ASUSY can be in agreement with LEP data. Ref. [H] also 
studied the issue of the Naturalness of the model, finding that this kind of models typically 
do not need to be tuned. The absence of tuning is somehow expected because the major 
source of fine-tuning in the MSSM is due to the need to generate large loop corrections to 
raise the Higgs mass above mz- Including the tree-level contribution coming from the large 
coupling A, this need is no longer a concern and one is not forced to push the model to an 
unnatural region of its parameter space because of the LEP bound on the Higgs mass. In this 
sense ASUSY is a remarkable candidate for a natural supersymmetric theory of Electroweak 
Symmetry Breaking. 

Furthermore it has been shown that in ASUSY a singlino-like lightest supersymmetric 
particle (LSP) can be a weakly interacting massive particle (WIMP) dark matter candidate 
with the correct thermal relic abundance. Additionally, ASUSY has a strikingly different Higgs 
sector with respect to the one of the MSSM and the one of the perturbative NMSSM, which 
leads to testable distinctive signals for the LHC [15] . 



An additional virtue of the NMSSM, if compared to the MSSM, is the possibility of generat- 
ing the Higgsino mass term \x dynamically. Indeed if the singlet S takes a VEV, the interaction 



in eq. (1.3) generates an effective mass for the Higgsinos jjl = X(S). In the NMSSM the VEV 
of the singlet depends on the same soft masses that trigger the VEV of the doublets. There- 
fore the \x term is generated by the same dynamics that breaks the EW symmetry, rendering 
evident why jjl is of the same order of the EW scale instead of being zero or of the order of 
some other energy scale that characterizes the (unspecified) UV theory. 

In this sense the solution of the /i problem in the NMSSM is very economical and consti- 
tutes a testable alternative to other mechanisms that generate the \x term through the same 
mechanism that generates the soft masses [HJJ [HI [181 HH] [2Q] [j 

Previous works on ASUSY [14] did not try to study the dynamical generation of \x and 
simply put the \i term by hand in the superpotential. In this work we shall address the issue 
of the dynamical generation of the /x term in ASUSY, which constitutes an important piece of 
information to complete the current picture of supersymmetry without a light Higgs boson. 

In our work we will discuss analytically the relation between the mass term jjl generated 
thanks to the interaction \SH\ ■ H2 and the other mass scales of the model. In particular we 
will discuss to what extent the dynamical generation of the \x term through the VEV of the 
singlet scalar requires the doublets H\ and H 2 to be mixed with the singlet. Because of the 
necessity of non-negligible mixing between doublets and singlet, we will generalize the analysis 
of [H] including the effects of the mixing where appropriate, as for instance in the EWPT 
analysis, in the consideration of limits from the direct searches of dark matter, and in the 
phenomenology at the LHC. 

A recent study of the NMSSM in the large A regime has been performed in Ref. [22], where 
the importance of a dynamical generation of a \x term is not highlighted and a numeric scan 
of the parameter space allowed by the many experimental constraints is performed. Contrary 
to this work, we will pursue an analytic approach as much as we can. Moreover the set of 
constraints that we will consider will be slightly different with respect to |22J . Notably, we 
will not require the NMSSM to provide an explanation for the current discrepancy between 
the experimental value and the SM prediction of the g — 2 of the muon. Besides, we do 
not impose the thermal production of lightest neutralinos to account for the observed relic 
dark matter abundance, as we content ourself to not overdose the Universe with the lightest 
supersymmetric particle. At variance with [22] we will take into account limits coming from 
direct dark matter searches through the elastic scattering of a weakly-interacting massive 
particle on a nucleus. 

Our paper is organized as follows. In the first part, we present the model and its theoretical 
predictions concerning the spectrum, the generation of the \x term and the level of fine-tuning 
required. In particular in Section [2] we describe the model and fix our notation. In Sections [3] 
and [4] we discuss the conditions for the global stability of the potential, for the preservation 
of the CP invariance and for a realistic breaking of the electroweak symmetry. In Section 

1 As a matter of fact explicit models of supersymmetry breaking generically have difficulties [TB] to generate 
correctly the fi term and special solutions for the generation of // are needed. In this sense the NMSSM appears 
more suitable for an economic dynamical generation of /1. Indeed it has already been considered the possibility 
to generate \x from an NMSSM superpotential with soft masses generated by gauge mediated supersymmetry 
braking [21] • 



[5] we summarize the constraints on the scalar potential for a correct electroweak symmetry 
breaking and for the preservation of the CP invariance, showing which is the allowed subspace 
of parameter space. In Section [6] we discuss the spectrum of the model. In Section [7] we 
analyse the generation of an effective \i term. In Section [8] we discuss the fine-tuning of the 
model and the predictions on the sparticle masses coming from Naturalness arguments. 

In the second part of the paper we discuss the experimental bounds that the theory has 
to satisfy. In particular, in Section [9] we discuss the limits from LEP direct searches. Then, 



in Section 10 we study the indirect constraints coming from Electroweak Precision Tests. In 
Sections HT1 and H2l we discuss the relic abundance of the LSP of the model and its detection 
in current experiments. 



Finally in the last part of the paper (Section 13) we discuss the phenomenology of the 



model and its possible signatures at the LHC and in Section 14 we give our conclusions. 



2 The model 

A recent review on the NMSSM can be found in [23], but for sake of completeness and to fix 
our notation, we give here some details on the model. We start considering the most general 
theory with SU(3)c x SU(2) L x U(l)y gauge invariance that contains, in addition to the Higgs 
fields of the MSSM, a singlet chiral field S. The most general superpotential that we can assign 
to the theory is W gen ($ l ) = {iE x ■ H 2 + fS 2 + XSHx ■ H 2 + §S 3 , where H x and H 2 are the 
Higgs superfields with —1/2 and 1/2 hypercharge respectively. This superpotential generalizes 
the superpotential of the MSSM and still contains dimensionful parameters such as \i and M. 
Although these parameters are protected by holomorphy, their presence is typically seen as a 
problem, the "/x problem" , because their size has to be fixed by hand, and in particular the 
\x parameter should be fixed to a value closed to the Z boson mass. To address this issue we 
impose a continuous R-symmetry such to avoid all the terms with a dimensionful constant. 
Consequently the theory is determined by 

WWssm^) = \SH X ■ H 2 + ^S 3 . (2.5) 

This superpotential plus the corresponding soft supersymmetry breaking potential 

V soit = m\ \H X \ 2 + m\ \H 2 \ 2 + /4 \S\ 2 - (AXSH^ + G^-S 3 + h.c.) , (2.6) 

o 

and the D-terms Vn coming from gauge interactions fully define our model. 
The D-terms of our model are exactly the same of the MSSM: 

Vy = \gl{\H 2 \ 2 -\H l \ 2 y, (2.7) 

V 2 = \g\ {h\TE x + E\TE 2 \ 2 , (2.8) 



where T l = y. Using the well known identity for canonical generators of SU(2) Yli^ab^cd 
26 a dSbc — S a bS c d, we can write again the total gauge potential as 



V D = V Y + V 2 = l - {g\ + g 2 ) (|# 2 | 2 - {H,] 2 ) 2 + \g\\E{E 2 \ 2 . (2.9) 

The total scalar potential of the theory will also include several terms involving the squark 
and slepton fields, but for our purposes it will not be necessary to deal with them. In the 
following we will assume for simplicity that all the parameters of the Higgs-Higgsino potential 
are real. 

In addition, we conventionally choose the couplings A and k to be positive. 

3 Stability of the scalar potential 

The first requirement we have to impose to the scalar potential is that it is bounded from 
below. At large values of the fields the quartic part of the potential dominates and, since its 
coupling is not negative, the only way to destabilize the potential is to make the quartic part 
vanishing along some direction. However, the quartic potential of the NMSSM is the sum of 
D-terms and F-terms which are each positive definite and, as we will recall in the following, 
they cannot vanish at the same time, if not for the trivial field configuration. As such, the 
global stability of the potential is guaranteed by the supersymmetric structure of the theory 
and no constraint for the soft terms emerges. 

More in detail, the global stability of the potential can be shown as follows. The quartic 
part of the potential gets contributions from the F-terms and D-terms, V^ = Vp + Vp. Since 

2 2 

9$ — 12i \Fi\ ■> the F-term part of the potential vanishes only if all the Fi vanish. 
This condition is satisfied, in general, along the direction 



v F = J2 i 



5 = 0, #i-# 2 = 0, (3.10) 

and for k = also along the direction of generic S and 

fli = 0, # 2 = 0. (3.11) 



Vp is explicitly given in (2.9) and is a sum of positive terms. Requiring each term to 



vanish, one finds that Vp vanishes only along the non-trivial direction 

#j# 2 = 0, (3.12) 

|#i|-|# 2 | =0. 



Consequently, for k ^ the function V^ vanishes only if both the conditions (3.10) and 



(3.12) hold, and this is only possible for the trivial configuration 

H 1 = 0, #2 = 0, S = 0. (3.13) 

Namely, in the k ^ case there is no non-trivial direction in field space along which the whole 
quartic potential vanishes. Consequently the potential is always positive at large values of 
the fields, and there are no constraints on the parameter space resulting from the condition 
of stability of the potential. 



In the k = case, instead, there is the additional F-flat direction of (3.11). Along this 



direction the condition for vanishing Vp (eq. (3.12)) is always valid and therefore the whole 
quartic potential vanishes. As such, the large field behaviour of the potential along this 
direction is dictated by the soft terms. Requiring the potential to be positive for large field 
values yields the condition /z| > 0. 



4 The minimum of the potential 

The minimum of the potential has to be a stationary point, therefore the extremal point 
conditions with respect to Hi = (H®, H®, S) must hold in a non-trivial point (t>i,t>2, s) 







V 



q h , v \ H o= Vl< HO=v 2 ,S= 







(4.14) 



which are equivalent to 



cos 2 (3 
X 2 v 



2„2 



m 2 + m z/ 2 + ^ S 

m\ + m 2 . + m| + 2A 2 s 2 

,2 ,^2 



., ,, 2sX(A — ks) m 2 — m, 9 

""" _ H ~„ +™> 2 Z , 



and 



sin 2/3 cos 2/3 

4fc V - 2Gks 2 + 2s (v 2 X(k sin 2/3 + A) + /x|) - Av 2 X sin 2/3 = 0. 



(4.15) 
(4.16) 

(4.17) 



The three conditions (4.15)- (4. 17) will be used as relations to trade v, tan/3 and s for the 



soft parameters mi, mi and \is- 

For fc = 0we can write explicitly the solution of the last equation for the VEV s 

AXv 2 sin 2/3 



2 (/4 + v 2 X 2 



(4.18) 



Instead in the general case of k ^ there are three different solutions of (4.17). As in the 



following we will impose the CP invariance of the vacuum, we will take the only real solution 
for s. This, however, will not be reported in the text, because of its quite complicated and 
not transparent expression. 

This stationary point has to be a minimum and therefore we will require the Hessian to 
be positive definite. In addition we want it to be a global minimum, hence we require it to be 
deeper than the origin 

V(y 1 ,v 2 ,s)<V(0,0,0)=0. (4.19) 

This latter condition reads 



/4s 2 + k 2 s A - 2G^s 3 - ^v 4 sin 2 2/3 



2 2 

m z v _ 2 



cos 2 2(3 < 



(4.20) 



while the condition on the Hessian happens to have a rather complicated expression when 
k 7^ and therefore we do not give it explicitly here. 



For the potential to break the electroweak symmetry down to electromagnetism, we have 
to ensure that the neutral Higgs boson potential has a global minimum in a point with non 
vanishing VEVs v\ and v 2 - Additionally, to ensure electromagnetism is not broken we have to 
impose a vanishing VEV for the electrically charged Higgs boson. Furthermore we will require 
the absence of spontaneous breaking of the CP symmetry, which is a necessary condition 
for the theory to not have tachyons at the extremal point [21]. This requirement amounts to 
impose that the VEVs of the scalar fields at the minimum of the potential have zero imaginary 
parts and will allow to treat the real and imaginary part of the scalar fields as non-mixing 
fields such that their mass matrices will be respectively a 3-by-3 and a 2-by-2 matrices, as 
opposed to the general case where a 5-by-5 mass matrix is needed. 

The requirement of positive VEVs v\ and v 2 is given by 



, m| 



Xs(A - ks) > \I(A- ks)Xs tan /3 - X 2 v 2 sin 2 /3 f- cos 2£ (4.21) 



' Xs^^- - X 2 v 2 cos 2 (3 + ^ cos 2/3 
tan/3 2 



where s is expressed by the real solution of (4.17) 



The discussion of the conservation of the U(l) em symmetry is simplified in the basis where 
only one Higgs doublet gets a VEV, v. In addition we exploit gauge invariance to set to zero 
the charged component of this same Higgs doublet such that we are left with a single charged 
scalar field that we call <fr ± . The absence of a VEV for ^ is then expressed by the condition 



d 2 V 



30±30±t 
which yields the condition on the soft breaking parameters 

Xv 2 „ , m 



>0, (4.22) 

</>±=o 



A > — sin 2/3 + ks - ^ sin 2/3 . (4.23) 

2s 2Xs 

For the discussion of spontaneous CP breaking we take for simplicity all the parameter 



of the Higgs potential to be real, we assume the condition (4.23) for the conservation of the 



electromagnetism to hold, so that we can write the potential at the minimum as 

Keutrai = A a |£| 2 (|l??| 2 + |# 2 °| 2 ) + X 2 \H[H^ + m\\H^ + m 2 |tf 2 °| 2 + /4|S| 2 + A; 2 |S| 4 + 

7 9 9 

- (AXSH^H , + G-S 3 - XkS 2 H^H<? + h.c.) + ^±^ ((^l 2 - |iJ 2 °| 2 ) 2 . (4.24) 
3 8 

With a suitable U(l) rotation of the scalar fields we can cancel the imaginary part of H±, 
so that we can write the three Higgs fields at the minimum as 

ffi = («i,0), (4.25) 

S = se lB , 



where v\, v 2 and s are taken real and positive. 



Replacing these fields in the scalar potential (4.24), and requiring a minimum in corre 



spondence of the point = 9 = 0, we find that for k ^ the conditions 

3Gks 2 - 5k\ Vl v 2 s + 2A\v t v 2 > , (4.26) 

\k (Gs(A - ks) - 3A\ Vl v 2 ) > (4.27) 

must hold. In the particular case of k = 0, these conditions get simplified to the unique 
requirement of A > 0. 

5 Parameter space 

In this section we give explicit bounds on the parameter space arising from the the various 
constraints discussed in Section |4j Each point of the parameter space is fixed by coordinates 
A, k, tan/3, s, A and G that are in principle unconstrained, however we will focus on particular 
regions described and motivated in the following. 

As motivated in the Introduction, we are particularly interested in the large A regime that 
helps in pushing the Higgs sector beyond LEP direct searches. Ref. [2] shows that this can be 
done for A > 1 and a small tan/3. Remarkably this is precisely the regime where contributions 
to the precision observables from the Higgs sector are better fitting LEP data. In the same 
work, all the contributions to the T parameter were worked out in a particular model coming 



from a superpotential of the type in eq. ( 1.3 ) that contains the mass term for the singlet MS . 
This superpotential in the limit M — > gains a U(l) symmetry under which the two doublets 
transforms with the same charge and the singlet with twice the charge of the doublet. This 
symmetry is customarily called Peccei-Quinn U(1)pq symmetry and is also a symmetry of our 
superpotential in the limit k — > 0. In the case of Ref. [T4] the [/(l)pQ-symmetric limit M — > 
resulted in too large contributions to the EW oblique parameter T from the Higgsino sector. 
In close analogy, we will explore the breaking of the U(1)pq symmetry through the coupling k, 
with particular attention to the case of large k, where the symmetry is dramatically broken^] 
A representative point of the interesting region of the parameter space is then 

A = 2, tan/3 = 1.5, k = 1.2 , (5.28) 

which is singled out by the request of a moderate tan/3 and the maximal values of k and A, such 
that the theory will stay perturbative at least up to O(10 TeV) [2S]. Having fixed a region of 



the parameters space in the surroundings of the point (5.28), the remaining parameter space 
is spanned by the VEV s and the soft breaking parameters A and G. Since the soft breaking 
terms are naturally expected to be of the same order of magnitude, we decide to make the 
simplifying assumption 

G = A. (5.29) 

In fact, a large hierarchy between A and G is not possible since the requirement of absolute 



minimum given in (4.20) does not allow G ^$> A and the constraint from CP invariance of 



2 We remark that the large breaking of the U(1)pq guarantees that there are no nearly massless pseudo- 
scalar to worry about in cosmology. 
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Figure 1: Physical region of the /x, m#+ parameter space for parameters fixed according to 



eqs. (5.28) and (5.29). 



(4.27) does not allow G < i. Variations on the assumption (5.29) are conceivable but we 



checked that they do not lead to any important change in our analysis. 

Finally we decide to parametrize the two dimensional parameter space s, A in terms of the 
parameters [i and m#+ that are introduced through the relations 



/i = As 



m H+ 



m 



w 



AV + 



2(jl(A 



¥) 



sin 2/3 
For later convenience we also introduce the combination of parameters 



~ 2 



m 



2 I i \ Z 

m„+ - m w + A v 



2„.2 



°H+ 



W 



(5.30) 



(5-31) 



We will see in Section [6] that fi and m#+ are parameters with direct relevance for the spectrum 
of the model and indeed they also allow a straightforward representation of the bounds we 
have to impose. This is evident in Figure [l] where we show the physical parameter space for 
the particular choice of parameters in ( 5.28[ ) and (5.29). The lower-bound on m#+ is given by 
the condition of absence of spontaneous CP violation given in (4.27), that in the new variables, 
after imposing the condition G = A, reads simply 



77) 



H+ 



> 2X 2 v 2 + m 2 w 



(5.32) 



The other bound mainly shaping the triangle is given by the condition of absolute minimum 



of (4.20) that in the new variables has the approximate simple form 



m H+ < m w + 4- 



F 



+ 



rrir 



sin 2 2/3 tan 2 2/3 



+ 



(5.33) 



which we comment later in Section [71 



The union of the requirements (5.32) and (5.33) gives a minimal allowed value for /i on 



6 Spectrum 

The Higgs sector of the theory contains seven bosonic degrees of freedom. In particular, in 
the mass eigenstates basis we expect one charged Higgs, three neutral scalar fields and two 
neutral pseudo-scalars, that do not mix with the scalars, since we imposed CP invariance. In 
order to investigate the spectrum of these six particles, it is convenient to express the scalar 
potential with the Higgs fields expressed by 



S = s + 



S 1 + iS 2 

~7T~ 



Hi 



i al *\ ( V\ + 



g V2 «1 



hi \ 


i CT7r 2 


( ° 


V2 ) 


, H 2 = e V2 v 2 


U+vi 



(6.34) 



After some algebra, we obtain the mass matrices given in the following, where we have 



used sp = sin/3, cp = cos/3, tp = tan/3, s 2 p = sin 2/3 and // and rriA as defined in (5.30) and 



(5.31) respectively. For the scalar Higgs bosons, in the basis (hi, h 2 , Si) we get 



Ml 






■ 2 m 2 z + s 2 ^m 2 A cpSp (2v 2 X 2 — m 2 z — m 2 A ) /j,v(2Xcp + spk) — sicpv\ r ^ \ 



m 2 z sl + c\m\ 



:k 2 ,,2 
= A 2 



4^/i 2 - Gfn + AVs 



+ v/j(2\sp + kcp) 

2„,2 2 m 2 A _i_ Afc„2 



2/3 4aj 2 



U 



2 - S W ) 



(6.35) 



For the pseudo-scalar we use the basis (tt^ 3 \ S 2 ), where n^ is defined by tc^ = sin/37r}' 



(3) 



cos /3ir 2 and we get 



-7,2 



M 



in 



v\cpsprh 2 A 



3kvfi 



PS 



X 2 v 2 c 2 s 2 ^ - 3kXcp S pv 2 + 3Gfi 



(6.36) 



The charged Higgs boson mass is just equal to the parameter m H + 



Z Z Z \ Z Z i 

m H ± = m H + = m w - X v + 



MA - f/j) 
sin 2/3 



(6.37) 



For the neutralino and chargino sector we assume that the gaugino mass parameters Mi , M 2 
are large. In this case the only light chargino is a pure Higgsino and has exactly the mass 



m. 



For the neutralinos we choose the basis defined by 
i = ^(Hi-H." " l 



V2^ 
and the mass matrix reads 

/ 

M N = 



N 2 = — [H l if,) . X, = S. 



(6.38) 



(6.39) 



^ 2 X(cp - sp) \ 

-V ~^X(cp + sp) 

\^X(cp-sp) ~^X(cp + sp) -2\n J 







(6.40) 
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Figure 2: Masses (in GeV) of the neutral scalars in the plane /i, m#+ for parameters fixed as 



in (5.28): (from the left to the right) the mass of the lightest CP-even (m Sl ), the masses of 

). In all the 



'm S23 ), and the masses of the two CP-odd scalars 



m 



0-1,2 



the two heavy CP-even 

panels the overlaid yellow area corresponds to m Sl > 114 GeV and the overlaid green area 

corresponds to the parameter space where the SU(2) breaking vacuum is stable and CP is 

conserved. 
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Figure 3: Higgsino masses (in GeV) as function of the chargino mass [i for fixed parameters as 



in (5.28). The dashed line corresponds to mz/2 which is taken as limit from LEP (eq. (9.79)) 
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From this mass matrix one can see that there is a massless Higgsino state if 

At> 2 A 2 sin2/3 



and that 



i=l,2,3 



/' 



<= 2 



k 



/r(|ilh c 2 X> 



X 2 



(6.41) 
(6.42) 



All the masses given in this section are plotted in Figures [2] and [3] in the plane /i, m#+ for 
the choice of parameters in (5.28). Few comments are in order. Firstly, we note that, due to 



the large value of A, the spectrum consists of relatively heavy Higgs bosons with a lightest 
CP-even mass of roughly (200-300) GeV and, due to the large value of k, there is no light 
state in the CP-odd sector, in fact the lightest CP-odd has mass of few hundreds GeV. This 
shows how this realization of the NMSSM is rather at odds with the widely studied case of 
small A. 

Another interesting point is the fact that the requirement of a mass for the lightest Higgs 
in accordance with the LEP bound tends to clash with the requirement of absolute minimum 
of (4.20). This implies that the model cannot have an arbitrarily large value of \x because 



this would yield a negative mass squared for the lightest Higgs boson. The existence of a 



maximal allowed fj, can be understood taking the the CP-even mass matrix (6.35) in the large 



\x and large A limit and observing that all the diagonal submatrices have negative eigenvalues 
independently of the other parameters. This fact is very welcome in view of the need to 
generate a /i term of the order of rriz and we will study the consequences of this fact in 
Section [71 

For later convenience, we conclude this section fixing some notation. We call the scalar 
and pseudo-scalar mass eigenstates (from the lightest to the heaviest) Si : s 2: Sz and ai,02 
respectively and we will denote the Higgsino mass eigenstates as Xi?X2,X3 (still from the 
lightest to the heaviest). 

We introduce the combination of doublet fields that take VEV, h = hi cos (3 + h 2 sin /3, and 
its orthogonal one H = hi sin/3 — h 2 cos/3. The rotation U from the mass basis to the basis 
(H, h, Si) is defined by 

s 2 j = U I h J . (6.43) 

Analogously for the pseudo-scalars the rotation matrix P is defined such that 



a i 
a 2 



P 



TV 



(3) 



S-2 



(6.44) 



For the Higgsinos we call V the rotation matrix from the basis (Ni, N 2 , N 3 ) to the mass 



eigenstates Xm such that 

Finally, we introduce the rotation matrix R^ for a generic angle x 



(6.45) 
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(cos x — sin x \ 
sinx cosx ] , (6.46) 

l) 

such that the states N n can be written as 

N n = RT >S + £ R nr' 4) H, . (6.47) 

i=l,2 

7 Generation of the fi term 

One of the original motivations of the NMSSM is the possibility to dynamically generate an 
effective ft term in the Lagrangian at the Fermi scale. This possibility is very interesting 
because in the MSSM such a term has to be fixed by hand to a value close to mz or the model 
would be badly unrealistic. Indeed such a term gives mass to the chargino and is crucial to 
have a correct EWSB. On the other hand in the NMSSM an effective /i-term is dynamically 
generated by EWSB through the VEV of the Higgs singlet resulting in 

fx = \s. (7.48) 

This dynamical generation of /i in connection to EWSB is particularly welcome also in view 



of fine-tuning problems. In facts the minimization equation (4.15) can be rewritten as 



, 1 2 m 2 — m 2 , tan 2 (3 
2 z tan 2 (3 - 1 



/i + -m z 2"— — - — = , (7.49) 



which shows that in a natural theory \i ~ mz- Generating the \x term via EWSB relates 
the size of \x to that of the soft terms and therefore solves at the same time the issue of the 
presence of the term and the issue of its size. 

In the traditional approach to the /z-problem within the NMSSM one requires that the 
coupling A has to stay perturbative up the GUT scale, which means A < 0.7 at the EW 
scale [12j [13] . As such, to obtain a value of /i compatible with current searches and with the 
requirement of EWSB, one has to go in the regime s > v. In particular, to have As ~ mz 



when A — > one expects s ^$> v , so that the minimization equation (4.17) has an approximate 

stable solution [26J 

1 



s^-(G+^G*-8nl). (7-50) 

This approximate solution generates larger values of \x as one takes larger A and therefore 
seems to signal that the theory will be fine-tuned when a too large A is considered. However 



this is not worrisome because the approximate solution (7.50) cannot be valid for any large 
value of A as the approximation s ^> v breaks down. 



Indeed one can consider the minimization equation (4.17) in the large A regime for generic 
s and find an approximate solution 

l_^4sin2^ 
S "2A;sin2/3 + A' U ' bij 
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that in the limit A — > oo gives 



/'c 



A lim s 



-A sin 2/3. 



The existence of such finite limit is not surprising as we have already noticed in Section 
[6] that (jl cannot be taken arbitrarily large because of the incompatibility of the requirement 



of absolute minimum (4.19) and the requirement of positive masses squared of the CP-even 



scalars. Indeed one can be more quantitative and show that the presence of such maximal 
value of /i is a generic feature of the model due to the largeness of A and that the value of the 
maximal /i is linked to the maximal mass of lightest Higgs boson. 

The actual maximal allowed value of \i can be estimated observing that in the region close 



to the boundary defined by the condition of absolute minimum given in (4.20) and for k < A 
the condition m% > can be approximated by 



m 2 H+ — rn^y + X 2 v 2 > /i 



3 — k/X 3 . ., .. 

*w + 2 vMk/x - l - 



(7.52) 



This is a lower-bound on m H + that, for values of \i large enough, gets incompatible with the 



condition (4.20) for the minimum of the potential. The value of \x where the two conditions 



get incompatible can be estimated taking (4.20) in the approximate form given by the upper- 
bound on m H + of (5.33), yielding a relatively simple condition in terms of p = k/X 



ji < 



vAsin2/3 3(p-4)p + ^/8(p - l)(5p - 7) + 9 _ 3 
2 (p-5)(p-l) ^2 



vX sin 2/3 + 



(7.53) 



from which we observe that for generic tan/3 and generic k < X the model has a maximal 
allowed value of \x of the order of Xv. 

Furthermore, from the requirements of absolute minimum and of the absence of sponta- 



neous CP breaking given in (4.20) and (4.27), we can see that there is a minimum allowed 



value for \i that can be estimated taking the approximate condition in (5.33) in place of (4.20) 
giving 



/i 2 > 



2„,2 



X 2 v 



sin 2 2/3 



rrir 



cos 2 2/3 . 



(7.54) 



Taking together this result and the condition (7.53) for the absolute minimum, we find 



that the model is consistent only for values of /x within an interval that, neglecting the term 
mz cos 2/3/2, reads 

^sin2/3< / i<^sin2/3. (7.55) 

In the large A regime this relation automatically ensures that the chargino with mass jjl is 



above the LEP bound (see later eq. (9.77)) and, at the same time, gives an upper-bound for /j, 



linked to the mass of the lightest Higgs boson. In this sense eq. (7.55) shows that, specializing 



the generic ASUSY superpotential in (1.3) to a superpotential without dimensionful parameter 



gives a model where the /i term is phenomenologically acceptable and is necessarily close to 
m z , thus solving the so-called "/^-problem" . It is important to stress that eq. (7.55) holds 
only for k < X and that this is generically the case in the large A regime of ASUSY. 
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Furthermore, we can use (7.55) to establish a relation between the mass of the chargino 

(7.56) 



and the mass of the lightest Higgs 

which is of course of phenomenological interest for collider searches. 



m x + ~ m si , 



8 Naturalness 

As mentioned in the Introduction, one of the motivations to consider this model is the attempt 
to address the Naturalness problem of the MSSM in the Higgs sector. Therefore, of particular 
importance is the study of the level of fine-tuning required, in order to satisfy the various 
constraints on the parameters of the model. 

But what do we mean with Naturalness constraints! Speaking in broad generality, the 
Fermi scale is a function of the several dimensionful parameters a,,- of the Lagrangian: v 2 = 
We require, for a small variation of the parameters, also the variation of v 2 not to be 



large as well. In particular, for a fixed maximum amount of fine-tuning ^, we impose 



i 



A, 



a]dv 2 \ 



0; 



da 2 



< A 



(8.57) 



For the particular theory analysed in this paper, the set of dimensional parameters is given 
by a, = (/is, mi,m2, G, A). To compute the several logarithmic derivatives in (8.57), we have 
to consider again the conditions of minimization of the potential (4.15) and (4.17), which 



show the dependence of the Fermi scale (4.16) on the several dimensional parameters through 
tan/? and s, respectively. Therefore, we take (4.15) for the angle /3 and we replace then the 



value of (3 in eqs. (4.16) and (4.17), to eliminate their dependence on tan/3. Subsequently, we 



eliminate the dependence on the VEV v in the equation for s (4.17), using (4.16). From the 



such obtained eq. (4.17), it is then possible to compute the several derivatives of s with respect 



to the dimensional parameters of the Lagrangian ds/da 2 . Finally, for the level of fine-tuning 



A,, . we obtain 



A, 



a 2 dv 2 ( 



da 2 



dv 2 (ai,s) dv 2 (a i) s) ds 



da 2 



+ 



ds 



da 2 



(8.58) 



From our analysis, it turns out that, considering just the dimensional parameters mi, rri2, 
G and /j,s the level of fine-tuning required is in general very small in all the parameter space 
allowed by the conditions on the scalar potential analysed in Section [4] and by the LEP bound 
on the mass of the lightest Higgs given in (1.2). 



The Fermi scale is more sensitive to variations of the soft parameter A, but also the fine- 
tuning required on this parameter is still under control. For example, for the reference point 
in (5.28) it is always smaller than ~ 25 for a lightest Higgs in agreement with the LEP bound. 



Our result for the fine-tuning on A is shown in Figure |4j as a function of the mass of the 
lightest Higgs, for the representative point of the parameter space presented in (5.28), once 
that the mass of the chargino is fixed to two reference values: /i = 230 GeV on the left panel 
and /i = 400 GeV on the right panel. Indeed, fixing these two values for fi, we do not loose 
of generality, since in the region allowed by LEP the fine-tuning depends only mildly on the 
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Figure 4: Logarithmic derivative of the Fermi scale with respect to the dimensional parameter 



A as a function of the lightest Higgs boson mass m si , for the representative point in (5.28) 



and mass of the chargino fixed to 230 GeV (on the left), and 400 GeV (on the right). The 
blue area is the region allowed by all the constraints on the scalar potential (see Section \4 
and by the LEP bound on the lightest Higgs boson mass. 



value of ii. Interestingly enough, we notice that, increasing the mass of the lightest Higgs, the 
fine-tuning decreases considerably. Consequently, for both values of fx, Naturalness arguments 
favor the regions of parameter space with a heavy Higgs boson, even quite heavier than the 
LEP bound mass of 114 GeV. 

In general we can argue that, with large values of the coupling constant A, one can address 
the fine-tuning problem present in the MSSM Higgs sector and partially present also in the 
NMSSM with small coupling constants [28] . 

We could have already guessed this final result, just looking at the minimization conditions 



(4.15)-(4.17). In fact, considering the limit A — ¥ oo, the following relations hold 



cos/3 


-)■ 


Constant 


dv 2 
dA 2 


->■ 


2\ 2 A' 


dv 2 
ds 


-> 


A 
A' 


ds 


— !► 


1 



dA 2 2\A 



(8.59) 
(8.60) 

(8.61) 
(8.62) 



Inserting these limits in the equation for the fine-tuning A^ (8.58), it is obvious that the 



Naturalness conditions are easily satisfied in the limit of large quartic coupling A, indepen- 
dently on the value of the coupling " A 



i A more careful analysis of the large A regime shows that the Naturalness conditions are even easier to 
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It is also interesting to understand if the theory continues to be natural, when we allow 



only the coupling k to be large. Considering the limit k — > oo, the relations (8.59)-(8.62) 
change in 



cos — >■ Constant , (8.63) 

dv 2 

— - -¥ Constant , (8.64) 

oA 

dv 2 

— -+kn, (8.65) 

ds 1 , 



dA 2 2kA 



Inserting these limits in (8.58), it follows that the fine-tuning generally does not decrease, 
increasing the value of the coupling k. Consequently, the only way to make the theory natural 
is to have a large coupling A regime. 

8.1 Naturalness bounds on sparticle masses 

Using Naturalness arguments, we want to set upper-bounds for the masses of sparticles. These 
bounds are particularly relevant to understand the expected size of the contributions of the 
sparticles to low-energy processes like the well studied flavour transitions and LEP precision 
data. Moreover with these bounds at hand one can estimate the timescale for the observation 
of such states at the LHC. 

For the case of the stops-sbottoms, in the hypothesis of diagonal squark mass matrices, 
the basic observation is that the soft mass rrig which enters in the formula for the masses of 
the physical stop and sbottom 



m JL \lnr () + m 2 + m|cos2/3 ( - --sm 2 6 w ) . (8.67) 



,2 , .../I ,._..-) .< / ^ I l 0^2; 



m b L = \l m Q + m b +^|cos2/3 [ — - + -sm z 9 w ) (8.68 

affects the soft Higgs mass iri2 through the one loop renormalization group equation (RGE) [29] 

!M._L A?K + <) + ..., (8.69) 

where the ellipsis stands for terms not dependent, in first approximation, on the soft squark 
masses. One can integrate these equations up to the messenger scale A mcss , obtaining, at the 
leading log, 



satisfy in the particular case of k large, even if the results for the fine-tuning A a . do not dramatically change, 
depending on the value of k. 
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Figure 5: Upper bound on the mass of the stop (left panel) and of the gluino (right panel) 
for h-mess ~ 100 TeV an allowed fine-tuning of 10%, computed for the representative point in 
(5.28) and for the chargino mass \i equal to 400 GeV. The blue area is the region allowed 
by all the constraints on the scalar potential (see Section El) and by the LEP bound on the 
lightest Higgs boson mass. 



Srrir, 



8tt : 



\2/ 2 , 2 



i -^mess 



.70) 



To give an estimation of the bound on the masses of stops and sbottoms, one can simply 



assume the equality of the soft masses m^ 



m 



tR- 



Imposing, then, the Naturalness condition 



(8.57), with respect to the parameter m*, and using the approximate expression 



m 



Q 



dv 2 



v 2 dm\ 



m 



Q 



dm 2 . 



dv 2 



dm 2 dm 2 . 



(8.71) 



one can find analytically the bound 



rriA < v 



471-2 



sin/3 



dv 2 
dm'i 



In 



A„ 



(8.72) 



1 TeV 



The result for the bound on the mass of the stop for a messenger scale of 100 TeV and an 
allowed fine-tuning of 10% is shown in the left panel of Figure [5] for the point of the parameter 
space presented in (5.28) and jjl fixed to 400 GeVrl From the Figure it is evident that the 



upper-bound on the stop mass mz increases when one increases the value of the mass of the 



lightest Higgs, and, for a Higgs mass satisfying the LEP bound, the mass of the stop can be 
~ 550GeV for an allowed fine-tuning of 10%. 



It is clear from eq. (8.70) that in the case of the sleptons and the first and second generation 



of squarks we would have a bound that is too loose to be useful because of the smaller Yukawa 



4 There are no relevant changes in the curve for the bound on the stop mass, for n fixed to be equal to the 
second reference value (230 GeV) discussed in the previous subsection. 



couplings. However, extending eq. (8.70) to include two-loops gauge effects the recent analysis 
of |8j shows that in our model Naturalness bounds on the sfermion masses are around 10 TeV. 
Such bounds are of no practical impact for our analysis. 

Additionally, it is of interest to know the Naturalness bound on the gluino mass rrig, which, 
contrary to the stops-sbottoms, contributes to 777.2 via a two-loop effect enhanced by the large 



QCD coupling constant a s 
details), one obtains 



Ml 

•17T 



Integrating the soft mass 7772 RGE (see Appendix 



B 



for 



5m\ 



a 2 A 

i/s \ 2 2 1 i 1 mess 

4^ A ^ ln Piw 



In 



A„ 



ITeV 



(8.73) 



Imposing the Naturalness condition (8.57), with respect to the parameter rrig, and using 
the approximate expression 



mj dv 2 



v 2 dm 2 : 



m\ dv 2 d 



mZ 



dm 2 , dm 2 , 
z 9 



(8.74) 



one can find analytically the bound 



m g < 2ii 2 vsm(3 



A 



9t 



dv' 2 



dm\ 1Xi 1 TeV 



(l — In Amcss_"\ 

I- 1 m 1 TeW 



(8.75) 



The result for the bound on the mass of the gluino for a messenger scale of 100 TeV and 
an allowed fine-tuning of 10% is shown in the right panel of Figure [5] for the parameter point 



presented in (5.28) and for \x fixed to be 400 GeV. From the Figure it is evident that the 



upper-bound on the gluino mass m g increases when one considers increasing values of the 
mass of the lightest Higgs, and, for a Higgs mass satisfying the LEP bound, the mass of the 
gluino can be ~ 1.3 TeV for an allowed fine-tuning of 10%. 

Our study of Naturalness yielded two interesting results. Firstly, the model is less tuned 
when the lightest Higgs boson is heavier. Secondly, the bounds on third generation squarks 



and gluinos are rather loose. In fact, we will see in Section 10 that LEP phenomenology is 
rather unaffected by the presence of the superparters of top and bottom. Similarly, the bounds 
on the masses of the sparticles involved in flavour transitions are quite loose and can naturally 
yield only negligible contributions with respect to the SM. 

These two points are in general quite interesting because they show how the model, in the 
region where the Higgs is quite massive, is able to settle the tension present in the MSSM for 
the masses of stops and gluinos. In fact, on the one side in the MSSM Naturalness arguments 
would require stops-sbottoms and gluinos to be light: if one does not want to allow a 1% or 
less of fine-tuning, stops, sbottoms and gluinos can not be at the TeV scale in the MSSM; 
on the other side, AF = 2 and AF = 1 flavour transitions would generally push to have 
heavier sparticles, since, otherwise the new physics effects in the flavour observables (b — > s'y 
for example) that are well measured and well in agreement with the SM prediction could be 
too large. On the contrary, in the model at study sparticles masses can be naturally raised 
at the TeV scale, and consequently the SUSY flavour problem can be addressed more easily 
than in the MSSM. Furthermore, the analysis of [S] confirms this result when first and second 
generation squarks are included in the discussion. 



19 



9 LEP direct searches 

Using the formulas for the masses of the Higgs bosons and Higgsinos given in Section |6j we 
can impose the bounds from LEP direct searches [30j EH U E2] 

m h > 114 GeV, (9.76) 

m x+ > 103 GeV, (9.77) 

m H ± > 79 GeV, (9.78) 

m xi > m z /2, (9.79) 

and exclude regions of the parameter space accordingly. For sake of simplicity, instead of taking 
the limits on the cross-section for neutralino production (summarized in |23J), we impose the 
tighter condition to have all the neutralinos with mass above mz/2. Similarly we are not taking 
into account the fact that the quoted limits apply to the MSSM and that in the NMSSM they 
could be less stringent. Thus the limits that we impose are rather conservative and serve the 
purpose of proving that there is a region parameter space that is certainly allowed by the 
current data. We are however aware that a larger region of parameter space could be allowed. 
The union of the constraints from vacuum stability found in Section [4] and the constraints 



from LEP is shown in Figure^ for the choice of parameters in (5.28) and two variations of 
A and tan/3. The strongest constraint on the parameter space comes from the bound on the 
neutralino mass given in (9.79[), which cuts a region in /x around the value given in (6.41) 



where the lightest neutralino is massless. On the other hand, we see that chargino, charged 
and neutral Higgs boson searches have basically no impact on the physical parameter space, 
once that the constraints on the scalar potential of Sections [4] and [6] are imposed. Indeed the 
maximal chargino mass that could be probed at LEP is less than the minimum allowed \x 



given in (7.54). Analogously, for the neutral and the charged Higgs bosons the typical mass 



in our model is significantly larger than the bounds and therefore only little restrictions arise 



from (9.76) and (9.78). 



10 Indirect constraints from EWPT 

It is well known that the NMSSM in the large A regime can have significant impact on the 
EWPT [2]. As such, we compute the new physics contributions to the oblique parameters T 
and S, coming from the new scalars and the new fermions of the model. In what follows, we 
perform the analysis of the EWPT in the S — T plane, with the experimental contours taken 
from [33] . 

Scalar contributions arise from the enlarged Higgs sector and from the sfermions. All the 
sleptons, the first and second generation of squarks are not restricted to be particularly light 
by Naturalness arguments (see Section [8]) and therefore their contribution can be neglected, 
assuming they are heavy. Third generation squarks are somehow special, indeed we have seen 
explicitly in Section [8] using Naturalness arguments, that they cannot be taken too heavy. 
Therefore the stop-sbottom system can potentially give a sizable contribution to the elec- 



troweak precision observables T and S which we compute in Section |10.2| The contribution 
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Figure 6: Six panels showing the area of the plane /i, m H + allowed by the constraints of a 
correct SU(2) symmetry breaking, CP invariance, and LEP direct searches. The left column 
is for A = 1.5 and the right column is for A = 2. In each column is shown, from the top to the 
bottom, the result for tan/3 equal to 1.5, 2 and 2.5. In all the cases k = 1.2 and for k ~ 1 the 
results are qualitatively unchanged. A solid blue line indicates the points where 2m Xl = m ai . 
The vertical cut of the green region is due to the LEP bound on m Xl . 
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arising from the enlarged Higgs sector with a heavy spectrum will be discussed in Section 10.3 
In the following section, instead, we discuss the contribution from the fermions that, in the 
limit of heavy gauginos, reduces just to the contribution of the Higgsinos. 

10.1 Higgsinos 

In the basis iV lj2i 3 the interaction Lagrangian of the Higgsinos reads 

£ = "f K*- (J^i - 7 M 7 5 ^ 2 ) + h.c. (10.80) 

+ f K (X- -fXT + Nrf-fN*) + |5 M (r Yx~ ~ iV l7 M 7 5 iV 2 ) . 

Therefore the contributions to T and S can be given using the loop functions F and A given 
in Appendix [A] by the expressions 

T= J2 E V^A^,m Xm )+ (10.81) 

i=l,2 m=l,2,3 



E E^ lm ^ 2n + v mV2 m ) 2 A(m Xm , -m Xn ) + 

m=l,2,3 n>m 

2 E (^imVam + Vamy^^m^j-mxJ, 



m=l,2,3 



5=^ ^(yi m y2n + yiny m ) 2 F(m Xm ,-m x J+ (10.82) 

m=l,2,3 n>m 

~ E (^lmVam + Vimy2m) 2 ^(m Xro , ~m Xm ) - F(/i, /i) , 



2 

m=l,2,3 



where V is the rotation matrix for the Higgsinos defined in (6.45). 

In Figure [7] we present the result as function of \x and the coupling k for several represen- 
tative values of tan/3 and A. The Figure shows that the generic value of the contribution to S 
is sizable but of the same order of the experimental uncertainty; on the other hand the value 
of T is more problematic and deserves some more discussion. 

The result of Figure [7] shows that T prefers small values of both A and tan/3, as expected 
from the dependence of the mass splittings of the Higgsinos on these two parameters. Moreover 
it shows that values of k close to zero or close the perturbative bound of ~ 1.2 give the 
best results, with the latter generically giving a better result than the former. This can be 



understood noting from (6.42) that the mass scale of the Higgsinos for fixed A goes like the 



product k/j, thus, in general, larger values of k and // tend to give a smaller contribution to 



T. However, for any finite value of /a, there is a non- vanishing value of k given by (6.41) 
that renders massless the lightest neutralino. In the regions of the li, k plane where (6.41) 
is satisfied, and of course in the vicinity of them, the contribution to T is enhanced by the 
presence of the light state. From this discussion it is clear that for k large enough the critical 
value of Li is pushed to be smaller than the minimal phenomenologically interesting ii ~ 
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away from the line where (6.41) is satisfied, a larger k gives a smaller 



100 GeV and that, 
contribution to T. 

This preference of the EWPT for large values of k gives further motivation to consider the 
regime of the NMSSM with large coupling k (see eq. (5.28)), namely the regime where the 
PQ symmetry is broken by a large coupling and hence all the states in the CP-odd sector are 
heavy. 



10.2 Stop and sbottom squarks 

We compute the contribution to S and T in the limit of diagonal squark mass matrices so that 
the interaction eigenstates for squarks coincide with mass eigenstates and the contribution to 
S and T are just 



T = 6A(m h ,m- bL ), 

S = F{m~ bL ,m~ bL ) - F{m- tL ,m lL ) , 



(10.83) 
(10.84) 



where the loop functions A and F are reported in the Appendix [A] and the masses of the third 



generation squarks are given by (8.67) and (8.68 ) respectively. In order to estimate the minimal 



effect of the stop-sbottom, we assume that the soft mass ttIq saturates the upper-bound due 
to fine-tuning constraints given in (8.72). 



Both the contributions to T and S are rather small compared to the experimental uncer- 
tainties. In Figure [8] we show the value of T in the plane fi, rriH+, for an allowed fine-tuning 
of 10% and the point of parameter space presented in (5.28). We find that, for a 10% fine- 



tuning and for a lightest Higgs boson in accordance with the LEP bound (yellow region in 
the Figure [8]), the contribution to T is always smaller than ~ 0.04, which is small compared 
to the contribution coming from Higgsinos. This feature is due to the relatively heavy stops 
and sbottoms, allowed by Naturalness constraints, and does not change significantly for larger 
values of tan/3. 



10.3 CP-odd and CP-even Higgs bosons 

For the Higgs sector we compute the contribution to T and S coming from the three CP-even 
states and the two CP-odd states. We compute the total value of S and T, taking as zero the 
value of the LEP Electroweak Working Group [31] minus the contribution that arises in the 
SM from a Higgs boson with a reference mass just above the LEP bound rrih SM = 115 GeV. 
This latter contribution is given by 



T ( m h SM ) = 3(A(m hsM ,m z ) - A(m hsM ,m w )) , 
s ( m h SM ) = F ( m h SM , mz) + m 2 z G(m hsM ,m z ) , 



(10.85) 
(10.86) 



with the loop functions A, F and G given in the Appendix |A~{ 

Subsequently, we add the contributions coming from the enlarged Higgs sector of the theory 
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Figure 7: Contribution of the Higgsinos to T and 5 in the parameter space /x, fe for different 
values of tan/3 and A. The left column is for A = 1.5, the right column is for A = 2. From 
the top to the bottom of each column are shown the results for tan/3 equal to 1.5, 2 and 2.5. 
Solid red lines with squared labels are the contributions to T, dashed black lines with round 
labels are the contributions to S. 
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Figure 8: Minimal contribution to T coming from the stop-sbottom sector in the plane /i, m#+ 



for the reference point of parameter space presented in (5.28), in the approximation of no 



mixing between ti, and £r, and with the stop mass saturating the upper-bound (8.72). The 
overlaid yellow area corresponds to m Sl > 114 GeV. 



T = Y,Vl T(m Si ) +Y t U? 1 A(m H+ ,m«) + 



10.87) 



i=l 

2 



t=l 



3 2 



+ J2 PfiA(m H+ ,m ai ) -J2J2 UlPlA(m sv m ai ) , 

3=1 i=l 3=1 

3 3 2 

5 = X)^(m„) + ^^^i^K^mJ - F(m H+ ,m H+ ] 

i=l j=l 3=1 



(10.88) 



where m Si and m ai are the masses of the scalar and pseudo-scalar Higgs respectively, the 



rotation matrices U and P are defined in (6.43) and (6.44) respectively, and the loop functions 
F, T, S and A are given in the Appendix [Aj 

Once tan/3, A and k are fixed, the mass spectrum and the mixing matrices still depend 
on the parameters /i, m#+. Therefore, in place of the customary plot of the position of the 
model in the S-T plane, we show in Figure [9] the contribution to T and S in the /i, ms+ 
plane for representative values of tan/3 and A. The contributions are generically well within 
the experimental uncertainty and sub-dominant with respect to the contributions from the 
Higgsinos sector computed before. 



11 Relic abundance of neutralinos 

The relic abundance of a heavy particle that decouples as a non-relativistic specie can be 
computed solving its Boltzmann equation [521 ESI EZ] • The resulting relic abundance is given 
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Figure 9: Contribution of the scalar sector to S and T in the plane /x, m#+ for different values 
of tan/3 and A. The left column is for A = 1.5, the right column is for A = 2. From the top to 
the bottom of each column are shown the results for tan/3 equal to 1.5, 2, 2.5. In all the cases 
k — 1.2 and for k ~ 1 the results are qualitatively unchanged. Solid red lines with squared 
labels are the contributions to T, dashed black lines with round labels are the contributions to 
S. The blue area denotes the region of the plane where all scalars have positive mass squared. 
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by the approximate formula 



ni2 1.07 x 10 9 GeV" 1 

Ml 2 ~ —-, : 9 , r , (11.89) 

^M P (ax f . . + bx 2 f0 2 V ; 



where a and b are constants related to the annihilation cross-section of the LSP, g* = 86.25 for 
m b <C Tf < m w is the number of SM degrees of freedom at the time of freeze out, M P the 
Planck mass and x/. . the normalized freeze-out point temperature Xf. . = Tf, ,/mLSP, with 
rriLSP the mass of the LSP, that in our model is supposed to be the lightest Higgsino xi- The 
freeze-out point can be found solving numerically the equation 

1 0m7g x M P m LS p{(Tv)^xj^. 

log — , (11.90) 



•"f.o. \/9* 

with g x = 2. 

The only relevant quantity from particle physics is the thermally averaged annihilation 
cross-section av of the LSP which can be expanded around its non-relativistic limit as 

(av)=a + bx. (11.91) 

As far as the cases where misp < f^w are concerned, the only available contributions 



to (11.91) are those coming from the annihilation \X ~* ff- I 11 general this process can 
be mediated by a s-channel exchange of a Higgs boson or a Z boson, and by the t-channel 
exchange of a sfermion. Since in this paper we are interested in the case of heavy sfermions 



see Section 8.1), we will only investigate the s-channel contributions. The most important 



contributions to the s-channel can be understood in terms of symmetry arguments. Indeed the 



annihilation cross-section, as expanded in (11.91), corresponds to a partial wave expansion of 



the annihilation process. Using CP properties of the mediators and the chirality structure of 
the intervening interactions, one can find that the s-wave annihilation can be mediated only by 
CP-odd scalars, while the p-wave receives contributions from both a CP-even scalar and the 
Z boson [38, 39J. The s-wave part of the cross-section is typically suppressed by the smallness 
of the Yukawa couplings involved, with the only exception of the on-shell production of the 
CP-odd scalar. The line in the parameter space where this condition is fulfilled is shown in 
Figure [6] as a solid line. Away from this line we can just repeat the analysis of the previous 
investigations of Higgsino DM in the large A regime [14] . In fact, making the identification 
M — > — 2/c/i/A, the Higgsino sector of our model is equivalent to the one studied in [Hj and 
therefore the resulting relic abundance is the same in the two cases. 

Above the WW and ZZ threshold new channels open and the LSP can now annihilate via 
a s-channel exchange of a Z boson or a CP-even Higgs boson, and a t-channel exchange of a 
chargino. Once again by symmetry arguments the s-channels contribute only to the p-wave, 
while the t-channel contributes to the s-wave and the p-wave. Also in this case the resulting 
relic abundance coincides with that one obtained in [14] . 



In Figure 10 , the resulting VLh is given as a function of /i for different values of tan /3 for 



both A = 1.5 (on the left) and A = 2 (on the right). The features of the curves of Figure 10 
are mainly due to the dependence of the cross-section on the mass of the LSP. In fact, when 
m^/2 < ttilsp < ww, the only annihilation channel is that one into fermions that is mediated 

27 




z 

Q 




H [GeV] 



H [GeV] 



Figure 10: Relic abundance of neutralinos for A = 1.5 (left) and A = 2 (right). The horizontal 
blue band corresponds to the 3<r interval from the 7-years WMAP result |1QJ. Dotted blue 
lines are for tan/3 = 1.1, purple solid lines for tan/3 = 1.3, dashed red lines for tan/3 = 1.5 
and finally dot-dashed green lines for tan (3 = 2. The vertical purple shaded area corresponds 
to m LS p < m z /2 in the case tan/3 = 1.3. 



by an off-shell Z boson, which yields a cross-section behaved as a ~ l/m 2 LSP . This means that 
the relic abundance increases with m^sp until new channels open. When annihilations into 
WW and ZZ are available, the relic abundance has to decrease accordingly with the opening 
of the phase-space for the new modes. Altogether the curves have a rise-and-fall shape with a 
maximum corresponding to misp ~ rriz and a maximal value being determined by the mixing 
coefficient that affects the annihilation into fermions. 

Generically the annihilation cross-section of the Higgsino DM is too large and the relic 
abundance is then too low to account for the observed amount of DM, Vtti 1 ~ 0.11 jlQ]. In 
spite of this general trend there are regions of the parameter space where m Xl < rriz and 
tan/3 < 1.5 where the LSP has a large singlino component [T31UT1H2] and the relic abundance 
is reproduced. 

The necessity of a very large component of singlino in the LSP to reproduce the correct 
order of magnitude of the relic abundance of the DM is well know [T^l W\\ W2\ and is due to the 
fact that a thermal relic with electroweak couplings should have mass around one TeV rather 
than few hundreds of GeV to satisfy the requirement of a correct abundance. In this sense 
the need for a large component of singlino in the LSP should not be regarded as an additional 
tuning of the model, at least with respect to the usual situation of supersymmetric models in 
which the dark matter is a weakly interacting particle with mass close to the weak scale. 

In the regions of parameter space that reproduce the correct relic abundance for the DM 
the model can be tested by direct searches of DM through elastic scattering which we describe 
in the following section. 
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12 Direct detection of the dark matter 



Direct searches of DM particles stored in the halo of our galaxy have been performed and 
several bounds on the properties of the DM exist [l3j S3]- The quantity probed by these 
experiments is the local DM density times the cross-section of an elastic scattering between 
a DM particle of mass M and the nuclei of the experiment [HJ H6]. The local density is 
typically assumed to be 0.3 GeV/cm , so that bounds are given directly on the DM-nucleon 
cross-section. Therefore in the following we will compute the cross-section for the scattering 
of the LSP of the model on a proton. 

As we have already discussed, in our model the sfermions are heavy and the LSP is a 
mixture of Higgsinos and singlino. Therefore the only particle that can mediate a DM-nucleon 
scattering is a Higgs boson interacting via the Yukawa couplings dictated by the superpotential: 



-C 



yuk = -j={S x H x H 2 + hxSH 2 + h 2 H x ~S) + -^SiS 2 + h.c. . 



These interactions contributes to the effective operator 



12.92) 



which mediates the spin-independent elastic scattering of the lightest neutralino xi on a nu_ 
cleus N probed in |33l EH] . 

The elastic cross-section at zero transferred momentum can be written as 



°si(XiP -> XiP) 



where the matrix element is given by 



1 



16ir(m p + m LSP ) 2 



\M\ 2 , 



(12.93) 



^— ' mi 



9xxs 



2m 2 p 

V 



12.94) 



where 



9xxsn 



A J2 2V la V lb U^ + V2kV 2 13 U^ , 



(a,b,c) 



12.95) 



and the indices (a, b, c) run over all the ordered permutations of (1, 2, 3), the matrix V is the 
matrix that brings the Higgsinos H i, H 2 , S to the mass eigenstates Xi an d the matrix U brings 
the scalar CP-even interaction eigenstates h\, h 2 , Si to the m ass eigenstates s m . In particular 
the V matrix i s rela ted to the matrix V defined in (6.45) through V = V t R ( - n ^^ and U is 
related to U of (|6.43|) through U = UR {/3+w/2 \ 



The effect of heavy quarks in the nucleon is taken into account accordingly to [17] and 
incorporated in the values of F u and Fj, which we take from chiral perturbation theory respec- 
tively equal to 0.11 ± 0.02 and 0.44 ± 0.13 [48J or, from QCD on the lattice, 0.14 ± 0.02 and 
0.23 ±0.01 0H]. 

The quantity affecting the most the cross-section asi is the mass of the lightest Higgs 



boson which is the scale that suppresses the operator of the interaction (12.92). In fact, the 
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spin- independent cross section can be estimated as 

1 ml ml /200 GeV\ 4 

ctsi^- J^-HT^cm 2 , 12.96 

lo7r v z m\ \ m si J 

that is of the order of the sensitivity of current experiments. Therefore direct searches of 
WIMPs can significantly restrict the allowed parameter space where Qh 2 < 0.1 [JO]. From 
the previous section, we know that for 1.5 < A < 2 such interesting regions are those with 
tan/3 < 1.5 and m#+ and /i in the regions outlined in Figure [6] 

Given the relevance of the mass of the lightest Higgs for the cross-section, the latter has 
a significant dependence on the parameters of the scalar sector [i and m#+. However none 
of the two is of direct significance for the experiments that, on the contrary, probe m^sp- 
Thus we will show our result trading \i for m^sp and we will fix some representative values of 
mtH+ taken in the range that is allowed by all the constraints on the scalar potential analysed 
in Section [4] The resulting DM-nucleon cross-section is given in Figure 11 as a solid and a 



dashed thin line, obtained using hadronic matrix elements from [49] and [48] respectively and 
assuming that the LSP accounts for the entire dark matter in the Universe. However the 
actual rate of DM-nuclei scattering in the model is typically reduced by the scarce amount of 



thermal relic neutralinos, therefore in Figure 11 we give also thick lines which correspond to 



the prediction of the model taking this reduction into account. 



From the thin lines in Figure [TTJ we see that the predicted cross section is typically above 
the lower-bounds from experiments. This shows that in this model is difficult to reproduce 
the relic abundance of the dark matter without violating the experimental limits on WIMP 
scattering. 

On the other hand when one takes into account the actual abundance of LSP computed 



in Section 11 the limits from direct detection experiments are not so stringent. The region of 
low niH+ is well below the limits because of a low nucleon-DM cross-section and a low relic 
abundance, while for larger m#+ the regions corresponding to m Xl > 60 GeV are typically 
excluded. In both cases the region that is compatible with direct WIMP searches indicates a 
preference for small of /i ~ (200 — 300) GeV. 

It is noteworthy to mention that for m Xl ~ (50 — 60) GeV the prediction is compatible with 
the cross section favored by the fit [50] to the recent claims [43] of observation of a DM signal 
in direct detection experiments. Typically the regions compatible with this intriguing hint of 
dark matter detection have relatively large fi and m#+ , however, the statistical significance of 
these claims is still relatively low and thus in the following we will concentrate on the region 
with lower m H + that is safely compatible with the direct searches of DM. 

13 LHC phenomenology 

After the discussion of the restriction of the parameter space from past experiments, we discuss 
now how the LHC can probe our model. The spectrum of the model has already been studied 



in general in Section pi Here we specialize the point selected in (5.28) fixing the remaining 
parameters to 

H = 240 GeV , m H+ = 520 GeV . (13.97) 
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Figure 11: Three panels showing \og l0 (asi) for the prediction of the spin-independent DM- 
proton scattering cross-section fixing tan/3 = 1.5. On the left for A = 1.5 and m H + = 500 GeV, 
in the middle A = 2 and m#+ = 550 GeV, on the right A = 2 and m#+ = 700 GeV. The 
dashed and solid lines correspond to the prediction obtained taking the values of the hadronic 
matrix element from [48] and from [19] respectively. Thick lines show the prediction of the 
model once the actual relic abundance is taken into account. The shaded green and cyan areas 
are those excluded by Xenon jH] and CDMS [43] respectively. 
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Figure 12: The spectrum of the model for fixed parameters as in eqs. (5.28) and (13.97) 



Solid lines correspond to neutral fermions, dashed lines to CP-even and CP-odd scalars, dot- 
dashed lines to charged particles. The two thin dashed lines correspond to m^/2 (in green) 
and 114 GeV (in red). 
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The resulting spectrum is given in Figure 12 and shows several interesting features. We notice 



that there is a light neutralino as LSP, xi> a second lightest neutralino, X2, rather close to 
the chargino, x + > °f mass H, and a significantly heavier \z- in the scalar sector we find a 
relatively heavy lightest CP-even Higgs boson, s%, that is rather close to its maximal mass 



(1.4) and to the second CP-even state, S2- Both these scalar states have a mass of order Xv 
and from the relation (7.55) we deduce that their closeness to the fermionic states x + an d 



X2 is in fact generic in the entire region of parameter space allowed by all the constraints 
discussed previously. The pseudo-scalars always have a heavy state close to the TeV and a 
light state of roughly few hundreds of GeV. This large separation between the pseudo-scalars 
is the result of a significant level-repulsion effect arising from the large mixing between the 
two states. Because of this effect, it is rather typical to have a spectrum with 2m ai < m si . 



The precise values of the masses in the point defined by eqs. (5.28) and (13.97) are 



m Sl = 307 GeV , m S2 = 352 GeV , m S3 



580 GeV . 



m 



XI 



66 GeV , m X2 = 255 GeV , m X3 = 609 GeV , 



m ai = 99 GeV , m a2 = 910 GeV . 



(13.98) 

(13.99) 

(13.100) 



From the spectrum it is evident that the regime of the NMSSM that we are considering 
is radically different from the case with perturbative couplings. Indeed in the spectrum we 
have Higgs and Higgsino states with masses of several hundreds of GeV, with the lightest 
state typically heavier than the Z boson. As such, a detailed analysis of the couplings and 
branching fractions is needed. This is performed in the next sections. 

13.1 Production of the new states 

We want to focus on the production of the scalar particles as they provide a direct handle 
on A, the characteristic coupling of the model. Furthermore the production of scalars is well 
studied in the SM and can be easily understood in terms of ratios between a coupling of a 
scalar of our model and the corresponding coupling of the SM Higgs boson of equal mass. 
In such a way one inherits all the studies available in the literature about (differential and 
inclusive) QCD effects in Higgs production. The ratios directly affecting the cross-section of 
the partonic production process 

(13.101) 



99 



Si, (Xj 



are those of the couplings of up-type quarks to CP-even and CP-odd scalars denoted respec- 
tively by £usi and £ tta . and given by 






(sin [3U* 2 - cos PU*,] 
cos/3 2 PJ\, 



(13.102) 
(13.103) 



where the matrices U and P are defined in (6.43) and (6.44) respectively. 

The reduced couplings squared £ ttSi for the the CP-even Higgs bosons are shown in Figure 



13 in the relevant region of the plane //, m#+ after that the constraints on the scalar potential, 
the constraints from direct searches at LEP, and from direct dark matter searches are imposed. 
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Figure 13: (Left column) Reduced couplings squared of the up-type quarks with the CP-even 
scalars as defined in ( J13.1Q2 ). (Right column) Reduced couplings squared of the SM vectors 
with the CP-even scalars as defined in (13.105). 
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The state s\ is in general significantly coupled to the up-type quarks and it has, over the 
majority of the parameter space, a cross-section of at least 50% of the one of the SM Higgs 
boson with the same mass. The state S2, on the contrary, has at most 30% of the cross-section 
of the SM Higgs boson of equal mass. The converse of this slight decoupling of s 2 is the 
relatively large coupling of S3, which can have a cross-section up to 15% of the one of the SM 
Higgs boson of the same mass. This quantitative analysis shows that the scalar CP-even sector 
consists altogether of three states all significantly coupled to up-type quarks. In particular 
we find remarkable to have such a large coupling for the heaviest state, S3. In general such 
mixing scenario seems favourable for the discovery of states beyond the lightest one. 

For the pseudo-scalars the value of the couplings to the up-type quarks are mainly deter- 
mined by the fact that the two mass eigenstates are almost maximal admixtures of singlet 
and doublet interaction eigenstates. As such, the two states nearly equally share a coupling 
equivalent to a fraction cos (3 of the Yukawa of the SM Higgs boson. In particular, for the case 



of the point of parameter space (5.28) and (13.97), we obtain £ ttoi ~ 0.16 and £ tta2 ~ 0.14. 



For the partonic production processes 

qq —7- VVqq — > qqSi and qq — ¥ VVqq — > qqSi , (13.104) 

we define the reduced couplings squared 

£ VVSi = Ul. (13.105) 



In this case, as we observe in Figure 13, the heavy states s 2 and S3 are more decoupled, thus 



their discovery in the vector boson fusion processes (13.104) is considerably more difficult than 
for a SM Higgs boson of the same mass. 

Finally, the fermionic states Xi;X2,X3>X + are produced through their gauge couplings, 
as in the MSSM. Similarly, the production of H + will occur through its gauge and Yukawa 
interactions, as it happens in the MSSM. 

13.2 Decays 



The ordering of the masses shown in Figure 12 and dimensional considerations on the decay 



width allows for a rough determination of the relevant decay channels of each state. We discuss 
some of them beginning from the lowest lying states. 

The state Xi is the lightest particle with negative R-parity and therefore it is stable. As 
well known its production results in large amount of missing transverse momentum. 

The next lightest state is typically a\ that, due to kinematics, can decay only in SM fermion 
pairs. Therefore its main decay modes are 

a^bb, (13.106) 

ai-»rf. (13.107) 

The lightest CP-even state, s\, has the usual decays into SM vector bosons 

s x -*ZZ, s x ->WW, (13.108) 
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Figure 14: The total width (left) and the invisible width (right) of the lightest CP-even 
Higgs boson in the region of parameters space compatible with direct constraints from LEP, 
dark matter searches and the stability of the potential. 



and the non-SM decays 



S\ — y CL\CL\ , 

si — > a\Z , 
«i ->■ X1X1, 



s\ -> XlX2- 



(13.109) 
(13.110) 
(13.111) 
(13.112) 



However, due to the fact that m. 



m. 



/i _ nu X2J 
the majority of the interesting parameter space. 



the decay Si — > X2X1 is n °t available over 



The state s\ has a total width of several tens of GeV, as show in Figure 14 . The contribution 



from the non SM-like decays is sizable and therefore the decays into W W + , which would be 
dominant for a SM-like state, is always subdominat, as shown by the branching fraction into 



W W + shown in Figure 15 



The reduced rate of the resonant production of W~W + through Higgs states that results 
from the lessening of branching fraction and the Higgs production cross section is the first 
evidence of the non-SM nature of the Higgs boson that we expect to show up at the LHC. 
The reduced production rate of vectors motivates the search of the Higgs boson in other final 
states. As show in Figure 15 when kinematically accessible, the decay into pseudo-scalars 
Si — > O1O1 is dominant. The process 



gg — t- h — > aiai — > rrbb , 



(13.113) 
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Figure 15: The branching fraction of the lightest CP-even Higgs boson into W~W + (left) and 
into aiOi (right) in the region of parameters space compatible with direct constraints from 
LEP, dark matter searches and the stability of the potential. 



is potentially interesting to clarify the reason of the reduced production rate of resonant 
W + W~ pairs. However we observe that the decay into pseudo-scalar is usually dominant, 
especially where the SM-like decay into vectors is more suppressed. As such, we find that the 
Higgs boson could be observed earlier in the non-SM decay S\ — > O1O1. A search for the Higgs 



boson in the channel eq. (13.113) might well be in the reach of the first few fb l of luminosity 



of LHC at 7 TeV center of mass energy. 

Among the non-SM decay modes of the state s\ the decay into two LSPs is particularly 
interesting because it results in an invisible decay of the Higgs boson of up to 10 GeV, as 



shown in Figure 14 The presence of a sizable invisible width is a sign of the large coupling 
A and therefore represents an indication of ASUSY. The invisible width can be measured 
reconstructing the line-shape of the Higgs resonance using decay modes in clean final states 
as Si —¥ ZZ — > 4£ and comparing the total width of the resonance with the widths of the 
observable channels. The measurement becomes easier when the decay into pseudo-scalars is 
absent as the branching fraction of the decay into ZZ gets larger. As such the measurement 
of the invisible width helps to cover the portion of parameters space where the decay into 
pseudo-scalars is not available. 

The states s 2 and s 3 are typically not heavy enough to decay into pairs of s\, thus the 
results of previous studies contained in [15] are not applicable to this case. We have seen 



in Figure 13 that s 2 and S3 are relatively decoupled from SM vectors, therefore the relevant 
modes are 

s 2,3 - > Oi\Z and s 2j3 — > a%ai , (13.114) 
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zz 


WW 


tt 


XiXi 


X1X2 


a±Z 


a\a\ 


T [GeV] 


Si 


0.088 


0.196 





0.090 





0.059 


0.568 


30.3 


S2 


0.004 


0.008 


0.002 


0.179 


0.027 


0.001 


0.782 


33.6 


S3 


0.023 


0.047 


0.039 


0.461 


0.013 


0.165 


0.255 


48.2 



Table 1: The branching fractions and the total widths of the CP-even scalars computed with 



parameters fixed as in eqs. (5.28) and (13.97). 



and those into fermions. For the latter case we notice that the modes 

S2,3 ->■ X1X2 and s 2 ,3 -> ti (13.115) 

are always kinematically allowed. However, because of the mass dependence of the partial 



widths and of the largeness of A and k, we expect the decays involving a\ of (13.114) to be 
dominant. 

The branching fractions and the total widths of the CP-even scalars computed with pa- 
rameters fixed as in (5.28) and (13.97) are given in Table [I] 



Due to the heaviness of the charged Higgs boson H + , the only two-body decay mode 
available for the chargino is 

X + -> XiW + . 

Similarly, the fact that m Sl > m X2 implies that the only relevant modes for \2 are the two-body 
decays 

X2 -* XiZ and \2 ~> X1O1 • 

The situation for a 2 is more involved as it has all the following decay channels: a 2 — > Sjd\, 
a 2 — > Zsk, a 2 — > fsAifsM , a>2 — > XtXj- Among the available decay modes, more than one 
involves large couplings and/or large final state multiplicities. Hence, a detailed computation 
of the partial width is needed to determine which channel actually dominates. 



14 Conclusions 

In this work, pursuing as much as possible an analytic approach, we closely examined the 
low energy phenomenology of a supersymmetric model with a strong coupling phase in the 
Higgs sector at some scale ~ (10 — 100) TeV. The presence of a strong coupling phase allows 
large couplings in the effective superpotential at the weak scale and helps to increase the tree- 
level mass of the Higgs boson, thus relaxing the tension between direct LEP searches and the 
requirement of a little level of fine-tuning of the model. 



The model under study is defined by the superpotential (2.5), which is a restriction without 



dimensionful parameters of the more general form ( 1.3 ). Models of the more general form have 



been shown to be realized as low energy effective superpotentials arising from UV complete 
models that are compatible with the unification of gauge couplings [51] . In this work we have 



shown that in the model defined by (2.5) one can address the so called jjl problem generating 
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dynamically the /z-term for the mass of the Higgsinos with a value of order mz in the same 
way as it is realized in the perturbative NMSSM. This, however, is non-trivial because the 
model under study is in a completely different regime characterized by large couplings A and 
k. 

Our quantitative study of the generation of /i shows that in this model the size of the \i 
term is connected to the mass of the lightest Higgs boson such that 

1 3 

-^"V ^ V ;$ 2 m *i • (14.116) 

The existence of such relation is welcome and not surprising. In fact, the generation of /i 
is a EWSB effect and the characteristic scale of the scalar sector responsible for the EWSB is 
\v, that is essentially controlling the Higgs mass. Hence the model is in a regime where the 
/j, term is generated always in connection with the mass of a heavy lightest Higgs and thus it 
is natural to have a \i term not too small compared to LEP direct searches though still close 
to the weak scale. We also studied the fine-tuning of the model in particular in those regions 
where /i is phenomenologically acceptable, by mean of the sensitivity of the VEV of the Higgs 
doublets to variations of the fundamental parameters and we found that the model requires 
the least fine-tuning when the mass of the lightest Higgs is large. Altogether we find that 
making the requirement of a heavy lightest Higgs one also gets a phenomenologically viable 
solution of the \x problem with low fine-tuning. 

In the second part of the paper, we checked the model against a great deal of experi- 
mental constraints and we found that the mass of the lightest Higgs plays an essential role. 
Some regions of the parameter space are already excluded by current data, but most of the 
more "natural" parameter space have still to be explored. The experimental situation can be 
summarized as follows. 

As far as the spectrum is concerned, the possibility of having a large quartic term A in the 
tree-level potential results in the fact that the Higgs bound from LEP can be easily satisfied, 
in particular for low tan/3. The relation between the lightest Higgs and the chargino masses 



eq. (14.116) explains why the mass of the chargino is generically above direct searches. The 
lightest pseudo-scalar and the lightest neutralino still have significant dependence on the choice 
of k. A small k would result in a tiny breaking of the PQ symmetry and thus in a "split" 
spectrum with a very light state in the CP-odd sector and all the rest of the scalars much 
heavier. On the other hand, a large coupling k would give a more coherent picture with a 
spectrum made entirely of heavy states. We pursued this latter possibility obtaining pseudo- 
scalars that easily fulfill LEP bounds. On the contrary, in the large k regime LEP limits on 
the lightest neutralino impose more severe bounds to the parameter space. 

Great care has been devoted to the study of the EWPTs, with the result that the larger 
contributions to S and T come from the Higgsino/Higgs sector. The model can be in good 
agreement with precision data and has a preference for large values of k and small values of 
tan j3 that give rise to a heavy lightest CP-odd and a heavy lightest CP-even Higgs, respectively, 
and altogether support once more the idea of a spectrum entirely made of heavy scalars. 

We also studied the thermal production of the lightest neutralino and we found that 
generically the abundance of relic neutralinos is too low to explain the whole dark matter 
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in the Universe. Notable exceptions are the regions of parameter space where m Xl < rriz and 
tan/3 < 1.5, that can give relic abundance within the interval given by WMAP. These regions 
can be further probed considering direct searches of WIMPs. It turns out that it is generically 
not possible to account for the dark matter in the Universe without violating WIMP searches 
limits. On the contrary when the relic abundance of neutralinos is less that the WMAP ob- 
servation one can be in agreement with WIMP searches when the mass of the lightest Higgs is 
above 200 GeV. This corresponds to cut away almost all the parameter space where A < 1.5 
and give further support to the idea of a heavy lightest Higgs boson and a large A. 



The experimental bounds are summarized in Figure 16 where the limits from direct searches 
are shown in the plane /i , m#+ . The Figure also shows the contribution of the Higgsinos to 
S and T, which can be taken as representative of the overall contribution of our model to the 
EWPTs. Altogether we find that for A ~ 1.5 the limits from WIMP searches and LEP rule 
out most of the parameter space. For larger values of A, WIMP searches allow a larger fraction 
of the parameter space, but limits from LEP still give significant bounds; in particular the 
limit on xi reduces the acceptable range of \i and the EWPTs require tan/3 < 1.5. 

We did not discuss the issues connected to flavour physics, however we expect them to not 
pose any problem. For example one can estimate the contribution to b — > s'j and realize that 
is rather generic to have small contributions beyond those of the SM. As consequence of our 
assumption of diagonal squark mass matrices, the only new contributions to b — > sj would be 
the loop of chargino-stop and the loop of charged Higgs and top. These two diagrams must 
cancel in the supersymmetric limit [52] but even barring the possibility of a cancellation we 
expect the arising contributions to be small. Indeed, the particles involved in the loops are 
in general relatively heavy (see the Naturalness limits on those masses in Section [8]), and, in 
addition, the value of tan/3 suitable for our model is quite small. Indeed one can estimate this 
contribution taking the limit of large charged Higgs mass such that only the chargino loop 
is left. Given the small values of tan/3 that we consider the resulting contribution is always 
small compared to the SM loop of W and top. A recent analysis (53] shows that for the values 
of /i ~ 200 GeV of interest for us is enough to take the lightest stop mass above 400 GeV to 
be in good agreement with the SM prediction even in the large tan ,9 case. As discussed in 
Section [8] such values of the mass of the stop are perfectly natural in our model. 

Finally, we want to point out the future possibilities to probe experimentally our model. 
Firstly, future experiments on WIMP direct detection have the potential to probe largely the 
parameter space of our model. Furthermore new experiments will clarify the nature of current 
low statistical significance claims of the observation of WIMP scattering out of a handful of 
events. In view of further clarification from forthcoming experiments we limit ourselves to 
report that our model has regions of the parameter space where the mass of the LSP and its 
cross section with the proton do agree with the rates and recoil energies of these claims. 

Secondly, the other main candidate to probe the model is the LHC. Naturalness arguments 
allowed us to estimate that scalar top partners and gluinos might well be in the 1-2 TeV range 
and are therefore in the mid-term reach of LHC. There are no constraints which forbid gluinos 
and stops as light as the reach at the 2010-2011 run of LHC, but there is no preference for 
such lightness. This is the counterpart of raising the Higgs mass at the tree-level. 

The scalar sector has light states that might well be in the reach of early LHC runs. The 
lightest CP-even state is substantially different from a SM Higgs boson. In a large fraction 
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of the parameters space it decays into WW and ZZ only subdominatly, while it decays 
dominantly in a pair of pseudo-scalars a%ai, when kinematically possible, or to a pair of LSPs. 
It seems optimistic to reach a discovery in the standard channels into WW and ZZ within 
the scheduled 2010-2011 run. Indeed the reduced production rate of resonant W~W + pairs is 
the first signal of the non-standard nature of the Higgs boson that we expect to show up in 
standard searches. Rescaling the cross-section from Ref. [51] for the LHC at 7 TeV center of 
mass energy, we estimate that at least a hint of the process 

pp — > si — > aiai , (14.117) 

that results in final states of the type bbrf, 4r,46 [55] might be visible with few 1/fb of 
integrated luminosity at the 2010-2011 run. 

Furthermore, a sizable invisible decay width of the Higgs boson is a further signal of the 
presence of non-SM decay modes, as for instance the decay into two LSPs. 

It is also rather generic to have a significant production of the light CP-odd with main 
decay mode into bb pairs. However, this signal seems difficult to observe because of the large 
QCD background of bb. As such, it seems more promising either to search for the light CP- 
odd boson in decay chains like those of the lightest CP-even Higgs, or to look for cleaner but 
suppressed decay modes like the one into rf. 

In the long-term run one can search for the remaining states of the CP-even sector, that 
can be all produced significantly at the LHC via gluon fusion and decay in detectable final 
states. Also the fermionic states of the neutralino and chargino sector are in the reach of a 
longer run of LHC. Their discovery and the determination of their mass is a crucial step to test 
the model, as the mass parameter \x can be determined directly from their spectrum. These 
states can be searched in final states with on-shell or off-shell gauge bosons and Higgs bosons 
and missing transverse energy. The long-term run of LHC should allow for the search of other 
heavy sates like the charged Higgs boson and the heavy CP-odd. 
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A One loop contributions to S and T 

We collect here the one loop functions used through the paper. 

For a boson loop with internal masses mi and m 2 and coupling the the gauge boson W^ 
given by iW^I d^fa, we have 
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Figure 16: Summary of the experimental constraints for A = 1.5 (a) and A = 2 (b). In both 
the panels (a) and (b) the lower row shows the constraints on the plane /i, m#+ for fixed 
k = 1.2 and tan/3 equal to 1.5, 2 and 2.5 from left to right. The green area corresponds to 
all the LEP constraints and all the stability constraints except that for the CP-even scalars, 
which is represented separately by the yellow area. The blue area corresponds to a DM- 
nucleon interaction rate below current bounds. In the upper row we show the contributions 
to S (dashed line) and T (solid line) coming from the Higgsino sector as a function of \x for 
the same choices of parameters. ,-, 



2a em v 2 A(mi,m 2 ) 
1 



1 



A-r 



F(m 1} m 2 ) 



16vr 2 

1 
96tt 2 



m 2 + m 2 . 



m\m\ 



In 



77li 



mt 



m% 



m% 



(A.118) 



In 



A 4 



+ 



4m 2 m 2 



m\m 2 



m 



m 



2N2 



+ 



+ 



mf + m 



o - 3m 2 m 2 (m 2 + m 2 .) m 2 
mf — m 2 .) 3 m^ 



(A.119) 



Additionally, also the diagram in which the gauge boson W^ and the Higgs boson <fi prop- 
agate in the loop contributes to the parameter S, giving as loop function 
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For a fermion loop with internal masses m\ and 7712 and a vector coupling W^ij^^, we 
have 
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Differently, for an axial coupling, the result can be obtained by letting m\ — » —mi- These 
expressions are valid for both Dirac and Majorana fermions, with an extra factor of 2 in the 
case of identical Majorana fermions. 

B Re normalization group running at two loops 

At the two loop level, the mass of the gluino rrig contributes to the renormalization group 
running of the Higgs soft mass m 2 . In this Appendix we collect the main equations and 



assumptions used in Section 8.1 to get eq. (8.73) for the running of m 2 



At two loops, the main dependence of m.2 on the gluino mass is given by 



dm 2 
~df 



57T" 



(A 2 (m| + m\ R ) + h 2 ) + ^9 2 s m- g (2X 2 t m- g - \\h t - X t hf) . (B.123) 



The second term is a purely two loop contribution directly dependent on the gluino mass, 
instead, in the first (one loop) term, we have dependence on the squark masses of third 
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generation and on the third generation trilinear term ht- At one loop, they are related to the 
mass of the gluino through the equations 

dm% drr?r o 

dt dt W y ' 9 ' l ' 

g 2 s \ t m~ g , (B.125) 



dh t 2 ,> 



dt 3vr 2 

that can be solved, using the leading log approximation. Assuming then that at the high scale 
ht is small, if compared to its running, then it is easy to prove that numerically the most 



relevant contributions in (B.123) are those not involving the trilinear term h t . 



Replacing then, inside the RGE for m-i the squark masses, obtained integrating the dif- 



ferential equation (B.124), one can finally prove eq. (8.73) for the running of m 2 due to the 
gluino mass. 

Addendum 

After that this paper was accepted for publication, Ref. [56] appeared and showed the possible 
appearance of absolute minima with a structure of the VEVs different than what considered 
in this paper. Generically the extremal points discussed in Ref. [56] can be deeper than 
the electroweak breaking minimum discussed in our text in significantly large regions of the 
parameter space of the NMSSM, however in our regime the potential evaluated at our extremal 
point is safely deeper than in most of the extremal points considered in [56J . 

The only exceptions are the extremal points with vanishing VEV of the Higgs doublets 
and non vanishing VEV for the singlet that can become deeper than our extremal point in 
sizable regions of the parameter space. Around these minima the electroweak symmetry is 
not broken and therefore the points of parameter space where these minima are the absolute 
minimum should be discarded 

Imposing on our parameter space that the minima discussed in Ref. [56] are not deeper 
than the one that correctly breaks the SU(2) symmetry, we observe that for k ~ 1 or larger 
this constraint only mildly affects the parameter space, cutting off only regions where the 
Higgs boson mass is not far from the LEP limit. According to the results of Section [8] and 



Section 12 the regions excluded by this constraint are the most fine-tuned and are typically 



in tension with the bounds from direct detection of the dark matter. 

After that this paper was accepted for publication a new analysis of our same superpoten- 
tial, but with k < 0, has been presented in Ref. [57J. They find that the constraint found in [56] 
is much less restrictive for k < than for k > 0, in particular there is acceptable parameter 
space for both small and large k. The conclusions of Ref. [57] on the expected phenomenology 
of a scale invariant ASUSY for k < substantially agree with our results obtained for k > 0. 



5 One could still keep these points if the SU(2) breaking vacuum is metastable with a life-time longer than 
the age of the Universe. 
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